FRW Cosmological model with Modified Chaplygin Gas and 
Dynamical System 



Nairwita MazumdeiQ 
Ritabrata Biswa^l 
Subenoy ChakrabortjH 

Department of Mathematics, Jadavpur University, Kolkata-700 032, India. 

Abstract 

O ■ 

^vq , The Friedmann-Robertson-Walker(FRW) model with dynamical Dark Energy(DE) in the form 

of modified Chaplygin gas(MCG) has been investigated. The evolution equations are reduced to 
an autonomous system on the two dimensional phase plane and it can be interpreted as the motion 
of the particle in an one dimensional potential. 
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£jy 1 Introduction 

The exciting observational evidences [TJ [2]in the last decade raise a challenge to the standard cosmology. To 
incorporate the present accelerating phase of the universe within the framework of Einstein gravity, one has to 
introduce a non-gravitating type of matter with a hugely negative pressure (of the order of its energy density) 
\ called dark energy (DE hereafter). For the mysterious DE, there are only very weak constrains on its form of 

an equation of state [UE]. The most common candidate for DE is the cosmological constant A, which can be 
considered as a perfect fluid with equation of state (EoS) p = —p, p = A. But it is discarded due to inconsistency 
\Q . in the value of A from type la supernovae(SNIa) observations in comparison with the value of A interpreted 

as vacuum energy (i.e., the Planck mass scale). So usually, investigations in DE are modelled as quintessence 
scalar field [4 or field with barotropic equation of state [5] . But the transition from a universe filled with matter 
to an exponentially expanding universe does not necessarily require the presence of a scalar field as the only 
alternative. Subsequently, attempts were made to use an exotic type of fluid-the so-called Chaplygin gas(CG) 
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•— and then it is extended as genera 1 " - i rr, '' r '' 1 — Li 
further generalised to modified CG(MCG) with EoS 0[7] 



having the EoS p — and then it is extended as generalised CG(GCG) p = — ■Jr, B > 0, < a < 1. It is 
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V = lP-^_ (!) 

with 7, P > and < a < 1. 

This EoS shows a radiation era (7 = |) at one extreme(when scale factor a(t) is vanishingly small) and a 
ACDM model at the other extreme(when 'a' is infinitely large). So at all stages it shows a mixture. Also at an 
intermediate stage the pressure vanishes and the matter content is equivalent to pure dust. 

In the present work, we formulate the dynamics of FRW cosmology with MCG as the matter in it and the 
evolution equations are shown to be represented in two dimensional autonomous system by suitable transforma- 
tion of variables. The nature of the critical points are analysed by evaluating the eigen values of the linearized 
Jacobi matrix. The paper is organized as follows : Basic equations for FRW cosmology with MCG are presented 
and dynamical system is formulated in section [2j Both finite and asymptotic critical points are analysed in 
section [3] The paper ends with a short discussion at the end, in the section |4j 
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2 Basic Equations for FRW cosmology with MCG 

The Friedmann equations which governe the dynamics are given by (assuming 8nG = 1 = c) 

a 2 k 



and 



and the energy conservation relation reads 
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p + 3-{p + p) = 
a 



From the above field equations the acceleration has the expression 

- = ~(p + 3p) 

a b 

Using the MCG EoS in the above conservation equation and integrating we have 

1 



p(a) 



(7 + 1) 



B 



a' 1 



(2) 
(3) 

(4) 
(5) 

(6) 



where p = 3 (7 + 1) (a + 1) and c is an integration constant. 
So the EoS parameter can be writen as 



u(a) = 



p 7c — Ba^ 



p c + Ba^ 

To form an autonomous system from the dynamical equations let us introduce a new variable 



(7) 



x = a 

then we have 

d = x 1 

(8) 

i = -f{l + 3w(o)}p(o) J 

Thus set of equations ([8]) form an autonomous system in the phase plane (a, x). Further using the phase space 
variables the field equation ((5J) can be written as 



-x 2 + V{a) 



with 



V(a) 



a 2 p(a) 



(9) 



(10) 



Equation © is termed as the first integral of the solution of the autonomous system. Also geometrically 
solutions of the dynamical system © should lie on the curve given by equation © on the phase plane. 



3 Critical points : Analysis of Dynamical System 



The critical points on the phase plane have coordinates 



(0, 0) and a c = 



c (1 + 3 7 ) 
2B 
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(11) 



2 



sssssssssssssss 

//////////////////// 

//////////////// 

1111111111111111111 
■ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

\\A\\\\\\\\\\\\\\\\\ 

X VWWWNWWWW 



0.90 0.95 1.00 1.05 1.10 




Fig. 1(a) Fig. 1(b) 

Fig. l(a)-l(b) represent the variation of a — x corresponding the value of the parameters a 
3.4; c := 0.7; 7 := 0.3 



1.3; B :-- 



and both lie on the a-axis. It is to be noted that the critical points correspond to static universe(s). Moreover, 
the first one corresponds to big-bang singularity and so we will concentrate mainly on the second critical point. 
As critical point lies on the curve (0 on the phase plane so we must have V{a) = — | i.e., k = +1 is the only 
possibility. Hence we can say that the critical points lie in the portion of the phase plane represented by the 
trajectories of the closed model. Moreover, the first integral corresponding to flat model, i.e., + V(a) = 
divides the phase plane into two regions namely %- + V(a) > or < which correspond to open and closed 

models respectively. Further the boundary of the strong energy condition is the straight line a = c ^g^ 
parallel to the x-axis and the second critical point lies on the line. The left of this line is the region in the phase 
space corresponds to decelerating phase (strong energy condition is satisfied) while the accelerating phase of 
evolution (strong energy condition is violated) is characterized by the right hand region of the above line [8] . 

Now we shall examine the character of the second critical point by evaluating the eigen values of the following 
linearization matrix A of the autonomous system : 
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(12) 



at (a—a c ,x—0) 



So Tr{A) = and det(A) = i £ [ap(a) {1 + 3oj(a)}] [acfi) = \a c p ac ^\ a= ^ with p{a c ) = (i^)^'" and 

duj_ _ g/£ca^j(l+7) 
da (C+Ba^Y " 

Then the eigen value problem, namely det [A — XI] — 0, i.e., A 2 — A x Tr(A) + det(A) = simplifies to 
A 2 + det(A) = 0. 
Thus we have 

(i)real eigen values of opposite sign if det(A) < 0. 

(13) 

(ii)purely imaginary conjugate eigenvalues if det(A) > 
In the first case the critical point is a saddle point while it is a centre in the other case. As in the present case 
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det(A) < so we have critical point which is saddle in nature and hence it is unstable in character [9] as shown 
in figures l(a),(b) for different choices of the parameters. 

To complete the study of the dynamical system [9] we now examine the critical points at infinity. To do 
this we make the following coordinate transformation in the phase plane to cover a circle S' at infinity : 



1 x 

(a, x) — > (p, q) : p = -, q = — ;p = 0, — oo < q < +oo 
a a 



Thus the autonomous system (J5J) transfers to 



(14) 



p = -pq 
q=-£{l + 3u}-q 2 



and the first integral becomes 



2p z V 



= -kp 2 



The critical points on the circle at infinity for the autonomous system (|15p are 

Pc = 



q c = ±V^f{T+M 



(15) 



(16) 



(17) 



which shows that critical points at infinity exists only in the portion of the phase plane where strong energy 
condition is violated(i.e., in the accelerating domain). It should be noted that for the autonomous system (|15|) 
there are other critical points in the finite domain which are already taken in to account. 

The linearization matrix which characterizes the nature of the critical points (|17|) is given by 



A : 



•U{p(1+3w)} 





2q c 



(18) 



So clearly, Tr{A) = — 3q c and let det(A) = 2q 2 and the eigen values are — q c and —2q c ,i.e., eigenvalues are real 
and of same sign. So both the critical points are node and they are asymptotically stable when q c > and are 
unstable for q c < 



4 Discussion 

In this work we study the FRW comology with MCG as the matter contained. The evolution equations are 
reduced into an autonomous dynamical system with a suitable change of variables. The dynamical system has 
a first integral and solutions to the system should lie on the curve represented by the first integral in the phase 
plane. Also the first integral can be interpreted as the energy conservation relation of a particle moving in an 
one-dimensional potential. Except the big bang singularity, there is only one critical point for the system which 
is saddle type and hence unstable in nature. Also there are two critical points on the circle at infinity and 
are node type asymptotically stable or unstable depending on the sign of the coordinate of the critical point. 
Finally, it should be noted that strong energy condition has a vital role for existence and type of critical points 
both at finite domain and at infinity. 
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